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A wide variety of coupled harmonic oscillators exist in nature. Coupling between different oscillators allows 
for the possibility of mutual energy transfer between them and the information-signal propagation. 
Low-energy input signals and their transport with negligible energy loss are the key technological factors in 
the design of information-signal processing devices. Here, utilizing the concept of coupled oscillators, we 
experimentally demonstrated a robust new mechanism for energy transfer between spatially separated 
dipolar-coupled magnetic disks - stimulated vortex gyration. Direct experimental evidence was obtained by 
a state-of-the-art experimental time-resolved soft X-ray microscopy probe. The rate of energy transfer from 
one disk to the other was deduced from the two normal modes' frequency splitting caused by dipolar 
interaction. This mechanism provides the advantages of tunable energy transfer rates, low-power input 
signals and negligible energy loss in the case of negligible intrinsic damping. Coupled vortex- state disks 
might be implemented in applications for information- signal processing. 

The magnetic vortex structure is very stable as the ground state in nanoscale magnetic elements 1,2 . It is 
characterized by an in-plane curling magnetization (chirality) around and an out-of-plane magnetization 
(polarization) in the central (core) region. In isolated disks, applied magnetic fields or spin currents induce 
vortex excitations, among which a translational mode 3 " 7 exists in which the vortex core rotates around its 
equilibrium position at a characteristic eigenfrequency (v 0 = (Dq/Iu) 8 typically ranging from several hundred 
MHz to ~ 1 GHz. This so-called oscillatory gyrotropic motion occurs because the gyroforce of the vortex core is 
in balance with the restoring force acting on it, both of which forces are due to the competition between the short- 
range exchange and long-range dipolar interactions 8,9 . The rotational sense of the vortex gyration is determined 
by the polarization p (Refs. 4, 8) and thus is controllable by thep reversal, where p = + 1 ( — 1) corresponds to the 
upward (downward) orientation of the magnetizations in the core. When the angular frequency of a driving force 
is close to cl> 0 , the vortex-core motion and its switching can resonantly be driven even with low power consump- 
tion 10 " 12 . Such a vortex oscillation in the linear regime functions analogously to a simple harmonic oscillator. 
Recently the vortex structure is attracting keen interest owing to its potential application to nano- oscillators that 
emit microwave signals 9,13 . 

Moreover, individual vortex oscillators in two physically separated magnetic disks can be coupled via their 
magnetostatic interaction 14 " 20 . A displaced vortex core and its motion in one disk generate dynamically rotating 
stray fields, consequently affecting the potential energy of the other disk, in which vortex gyrations can be 
stimulated resonantly, and thus affecting each other. The relative vortex- core displacements in both disks as 
well as the disk-to-disk interdistance modify the interaction strength. These behaviours are analogous to those of 
coupled harmonic oscillators such as coupled pendulums or capacitively- coupled inductor-capacitor resona- 
tors 21 . Under the free -relaxation condition of such coupled oscillators, the kinetic and potential energy of one 
oscillator can be transferred to the other 21 . 

In the present study, we demonstrated reliably controllable dipolar-coupled vortex-core oscillations by a direct 
experimental observation of a complete energy transfer and all of the collective normal modes. This mechanism 
is a robust means of tunable energy transfer and information- signal transport between physically separated 
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magnetic disks, and provides for the advantages of negligible-loss 
energy transfer and low-power input signals using the resonant 
vortex excitations. 

Results 

Experimental observation of vortex gyrations in dipolar-coupled 
vortex oscillators. Figure 1 shows the sample, which contains several 
two-disk pairs, each of which consists of two Permalloy (Py: 
Ni 80 Fe 2 o) disks of the same dimensions. Note that both disks are 
positioned on the x axis, which is referred to as the bonding axis. 
The dipolar interaction between the two disks along this axis breaks 
the radial symmetry of the potential well of isolated disks. Applying a 
pulse current of 90 ns width and 2.5 ns rise and fall time along a 
single-strip Cu electrode (along the y axis) placed on only one disk 
(here denoted as disk 1), we shifted the core position in disk 1 to 
— 165 nm in the +y direction from the center position (x iy y Y ) = 
(0,0), and then allowed it to relax after turning off the field generated 
by the current pulse. During this free relaxation, the dynamics of 
the vortex gyrations in both disks of each pair were observed 
simultaneously, using spatiotemporal-resolved full-field magnetic 
transmission soft X-ray microscopy (MTXM). The magnetic 
contrast was provided via X-ray magnetic circular dichroism 
(XMCD) at the Fe L 3 edge (—707 eV). Measurements were made 
on the basis of a stroboscopic pump-and-probe technique 18 ' 22 ' 23 (see 
Supplementary A for details). In our earlier work 18 , the microscope of 
70 ps temporal and 20 nm spatial resolution was used to resolve the 
phases and amplitudes of both vortex-core positions, for the purpose 
of a feasibility test of such a state-of-the-art experiment. Figure lb 
shows the two pairs studied, that is, d int /(2R) = 1.05 and 1.10, and 




Figure 1 | Sample geometry of pairs of two vortex-state disks and initial 
ground states, (a) Scanning electron microscopy image of sample 
containing several two-disk pairs. Each pair contains two Permalloy (Py: 
Ni 80 Fe 2 o) disks of the same dimensions, as indicated in the inset. The center- 
to-center distance normalized by diameter d int /(2R) varies, such that 
di nt /(2R) = 1.05, 1.1, 1.15, 1.2. Each pair is separated 5.6 /mi from the 
neighbouring pairs. The insets show schematic illustrations of the sample 
with the indicated dimensions and the employed current pulse of 90 ns 
width, 2.5 ns fall-and-rise time, and 7.0 X 10 6 A/ cm 2 current density, (b) Fe L 3 
edge XMCD images of both disks for two pairs, d int /(2R) = 1.05 and 1.10, the 
states are represented by perspective-view simulation results, as indicated. 



includes simulation -perspective images of the initial vortex ground 
states in both disks. For the other pairs having larger d int /(2R) values, 
the gyrations in disk 2 were difficult to observe, due to the small 
deviation of the vortex core from the center position, caused by the 
rapidly reduced interaction strength. 

Figure 2 shows the experimental data on the vortex gyrations 
in the two disks in both pairs: d int /(2R) = 1.05 and 1.10 in (a) 
and (b), respectively (see also Supplementary Movies. 1 and 2). 
Representative serial snapshot images of the XMCD contrasts are 
shown in each top panel. The trajectory curves of the motions of both 
vortex-core position vectors, X x = (x l5 y Y ) and X 2 = (x 2 , y 2 )> are 
plotted in the lower-right panels together with, in the lower-left 
panels, their oscillatory x and y components. Since a quasi-static local 
magnetic field (a pulse of sufficient length: here, 90 ns) was applied 
only to disk 1, the core position in that disk was shifted to (x l5 y x ) = 
(0, +165 nm)for[p 1 ,C 1 ] = [ + 1, +1] andtofe, = (0, -165 nm) 
for [pi, Ci] = [ + 1, —1]. The opposite sign of the y displacement 
in disk 1 between d int /(2R) = 1.05 and 1.10 is attributable to the 
opposite chirality of disk 1 in the two pairs. Once the pulse field 
was turned off (at t = 0, the reference), the vortex core in disk 1 
began to gyrate with decreasing oscillation amplitude starting from 
(*i> y\) = (0> + 165 nm), whereas the vortex core in disk 2 began to 
gyrate with increasing oscillation amplitude starting from (x 2 , y-i) — 
(0, 0). The vortex core in disk 1 gyrated counter-clockwise, indicating 
the upward core (p x = +1) orientation, whereas the vortex core 
in disk 2 gyrated clockwise, corresponding to the downward core 
(pi = —1) orientation (see Fig. 2a). Figure 2b shows similar char- 
acteristic oscillations of the given vortex state configuration. 
Micromagnetic simulation results (see Methods) are in excellent 
agreements with the experimental one, except for small discrepancies 
of their frequencies and absolute amplitudes. Experimental results 
show smaller oscillation amplitudes and frequencies, which might be 
attributed to the imperfection of the sample and weaker interaction 
between both disks due to a possible reduction in the saturation 
magnetization of the sample. 

Energy transfer between two dipolar-coupled vortex oscillators. 

The most important finding here is that the vortex- core gyration in 
disk 2 is stimulated by the vortex- core gyration in disk 1, not by 
externally applied pulse fields, but rather through the magneto - 
static interaction between the two disks. For the case of d int /(2R) = 
1.05, the decreasing amplitude in disk 1 started to increase again at 
the time t = 26 ns, whereas the increasing amplitude in disk 2 began 
to decrease (see Fig. 2a). For d int /(2R) = 1.10, the amplitude decay in 
disk 1 occurs at a longer time, t = — 45 ns (see Fig. 2b). Since the 
magnetic potential energy of a vortex depends on the displacement of 
the vortex core from the disk center 8 , a complete energy transfer 
between dipolar-coupled vortex-state disks can be observed 
directly, via the variations of the displacements of both vortex 
cores, that is, |X X | and |X 2 |> as shown in Fig. 3. Such a crossover in 
the vortex-core gyration modulation envelope between disks 1 and 2 
confirms that the energy stored via the displaced core by an input- 
signal pulse field is transferred from disk 1 to disk 2 through the 
robust mechanism of stimulated vortex gyration in dipolar- coupled 
vortex oscillators. The high efficiency of this energy transfer 
mechanism is evidenced by the low (almost zero) amplitudes at the 
nodes of the beating pattern of the coupled vortex oscillators. 
Note that the mismatch between the nodes of the vortex-gyration 
amplitude in one disk and the anti-nodes in the other disk 
results from the intrinsic damping of a given material (see 
Supplementary B). 

The total energy of undamped coupled oscillators, which is the 
sum of the first and second oscillator energies and the interaction 
energy between the oscillators, is constant. Magnetic vortices 
typically reside in potentials, which give rise to corresponding eigen- 
frequencies up to — 1 GHz. However, in principle these frequencies 
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Figure 2 | MTXM observation of vortex-core gyrations in dipolar- coupled Py disks for d int /(2R) = 1.05 in (a) and 1.10 in (b), compared with the 
corresponding simulation results. In each of a and b, the first row, are serial snapshot XMCD images of temporal evolution of vortex-core gyrations, 
starting from t = 8.8 ns (©) and ending at t = 19.7 ns ((g)) for a), and starting from t = 12.5 ns (©) and ending at t = 24.2 ns (®) for b). In each of a) and 
b), the second (disk 1) and fourth (disk 2) row are experimental results of the x and y components of vortex-core positions (left) from center position 
(Xy y) = (0,0), and constructed vortex-core trajectories (right). The third and fifth row indicate the corresponding simulation results. 



can be modified virtually in any desired frequency range, which is an 
important and advantageous aspect of the application of coupled 
vortex oscillators. As in every coupled harmonic oscillator system, 
the original resonance frequency of the non-interacting oscillator is 
either enhanced or reduced by the coupling 21 . Generally a symmetric 
and an antisymmetric mode having a lower and a higher frequency 
compared with the original eigenfrequency in uncoupled disks, 
respectively, appear. The frequency splitting between the symmetric 
and antisymmetric modes have a direct correlation with the inter- 
action energy, and can be deduced from the beating pattern of 
coupled oscillations (see Supplementary B). It is clear that such a 
modulation envelope crossover between disks 1 and 2, as demon- 
strated in Figs. 2 and 3, is the result of the superposition of two 
equally weighted normal modes of different frequencies excited 
in a given dipolar- coupled oscillator system. We assume that all of 
the experimental results are within the linear regime. The beating 
frequency Aco of the modulation envelope functions in disk 1 and 
disk 2 can be obtained by fitting the data of |X X | and |X 2 | to two 
orthogonal sinusoidal functions with an attenuation term, |X X | = 
|Acos(Acof/2)|exp(-/ft) and |X 2 | = |Asin(Acof/2)|exp(-/ft): 
Aw/2n = 20 MHz and 12 MHz for d int /(2R) = 1.05 and 1.10, 
respectively. These results agree well on the fact that this angular 
frequency splitting Aco takes place according to interaction between 
oscillators in a coupled system (see Supplementary C) and the mag- 
nitude of Aco varies with the interdistance between neighbouring 
oscillators. 

Another important parameter is which represents the energy 
loss during vortex-gyration-mediated signal propagation through 
the neighbouring disks. Based on Thiele's equation of vortex-core 
motion, the attenuation parameter is expressed analytically as 



P = —{DIG)coq for an isolated disk with the gyro vector constant G 
and the damping constant D (Ref. 24). This parameter can be rewrit- 
ten as P = \ql\2 + \n(R/R c )]coo with the intrinsic Gilbert damping 
constant a and the vortex-core radius R c = 0.68L ex (L/L ex ) 1/3 where 
the disk thickness L and the exchange length L ex = {lA ex / Mf) 1 ^ 2 
with the exchange stiffness A ex and the saturation magnetization M s . 
For the case of experimentally obtained value of a = 0.01 ± 0.002 
(Refs. 25-27) for a given Py material and the eigenfrequency 
co 0 = 2n X (157 ± 3) MHz for isolated Py disks of R = 1.2 /mi 
and L = 50 nm (see Methods), the value of ft = 2n X (4.81 ±1.05) 
MHz is found to be in good agreement with the values of p , In X 
(4.2 ± 0.7) MHz and In X (3.0 ± 0.7) MHz for d int /(2R) = 1.05 and 
1.10, respectively, which are obtained from fits to the experimental 
data shown in Fig. 3. Such good agreements are strong proof that 
energy loss through vortex-gyration-mediated signal transfer is 
determined by the intrinsic damping constant and dimensions of a 
given material. This means that further reduction of attenuation can 
be achieved by material engineering. For example, for a relatively low 
intrinsic damping constant, a = 0.0023 for NiMnSb, ft is reduced 
drastically to 2n X (0.8 ±0.1) MHz (see Supplementary D). Some 
Heusler alloys are reported to have extremely small values of a, for 
instance, a = 0.00006 for Co 2 MnSi as found from an ab initio cal- 
culation 28 and oc= 0.001 for Co 2 FeAl as obtained from ferromagnetic 
resonance measurement 29 . 

Normal modes representation of coupled vortex oscillations. To 

directly obtain the value of Aco from the experimentally observed 
vortex gyrations in the coupled system, we transformed the x and y 
component oscillations in each disk into two normal-mode oscillations 
on the basis of normal-mode coordinates (Nl x , Nl Y ) — (*i +x 2 ,yi —y-i) 
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Figure 3 | Comparison of vortex-core displacement variations versus 
time in both disks for d int /(2R) = 1.05 and 1.10. The open and closed 
circles represent the experimental results for disks 1 and 2, respectively. The 
thick solid curves are the results of the fits to I Xi I = I Acos( Acot/2) I exp( — fit) 
and IX 2 I = IAsin(AcL>*/2)lexp(-/fr), where amplitude A = 141±12 nm, 
attenuation parameter ft = 2tiX(4.2±0.7) MHz, Aw = 2tiX(19.3±2) MHz 
for diJ(2R) = 1.05 and A = 166±15 nm, ft = 2tiX(3±0.7) MHz, and 
Aw = 2n( 1 1 .2± 1 .5) MHz for d int /(2R) = 1 . 10. The vertical lines correspond 
to the rate of energy transfer between the two disks: z ex = \{2k/ Aw). 

and (N2 X , N2 Y ) = (xi—x 2 , yi+yi) for the case of pip 2 = —1. In the 
normal-mode representation, the experimentally observed vortex 
gyrations in real coordinates (x, y) were decoupled into two normal 
modes having corresponding single dominant frequencies, as shown in 
Fig. 4 (see also the two normal modes for a different value, d int /(2R) = 
1.10, in Supplementary E). Figures 4a reveals that the two normal 
modes were excited almost equally and that their amplitudes were 
damped monotonically, as in free-relaxation vortex- core gyrations in 
uncoupled disks. According to their FFT spectra in Fig. 4b, the lower 
and higher eigenfrequencies for the two normal modes were 145 ± 10 
and 165 ± 10 MHz, and thus, their frequency splitting was Aco/2n = 
20 MHz for d^HR = 1.05, which is in excellent agreement with the 
values obtained from the fit to the modulation envelope functions of 
\Xi\ and |X 2 | with the corresponding damping (Fig. 3). 

Moreover, we can decompose coupled vortex gyrations into the 
two normal modes of each disk with respect to the ordinary coordi- 
nates according to the relations of (x\, y\) m = |(Nlx, Nly), 
(xi>yi)m = \( m x, N2 Y ), (x 2 , y 2 ) m = \{N\ X , -N1 Y ), and 
(x 2 , yi) N2 = \( — N2 x , N2y). Figure 4c shows that for the higher- 
frequency Nl mode, the x positions of the vortex- core motions in 
both disks oscillate in-phase but their y positions do so out-of-phase, 
and vice-versa for the lower-frequency N2 mode. Since the vortex- 
core motion in one disk in the direction normal to the bonding axis 
(i.e., along the y axis) results in stronger stray fields (exciting vortex 
motion in the other disk) than it does along the bonding axis (the x 
axis), the dipolar interaction is dominated by the relative y position of 
vortex- core motions between the two disks. Thus, the interaction 
energy of the Nl mode having out-of-phase oscillation at the y 
positions is higher than that of the N2 mode having in-phase oscil- 
lation at the y position. Analogous to a general coupled harmonic 
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Figure 4 | Normal-mode representations of vortex-core gyrations in 
dipolar-coupled oscillators for d int /(2R) = 1.05. (a) Oscillatory core 
motions and (b) dominant frequency spectra. The two normal-mode 
coordinates are represented as (Nl x = Xi+x 2 , Nl Y =yi—y 2 ) and 
{N2 X = Xi~x 2 y N2 Y = 71+72) > corresponding to their trajectory curves and 
FFT spectra, c, The x and /positions of vortex-core oscillations in each disk 
for each normal mode (see the text). 

oscillator, the Nl mode with the antisymmetric y component oscil- 
lation works harder than the N2 mode with the symmetric one. As the 
result, the frequency of the Nl mode is higher than that of the N2 
mode, as noted earlier. More generally, in two dipolar- coupled disks 
with arbitrary p- and C-configuration states, the mode characterized 
by out-of-phase oscillation between C x y x and C 2 y 2 is higher- 
frequency antisymmetric, whereas that with in-phase oscillation is 
lower-frequency symmetric. This is owed to the fact that the inter- 
action energy of each mode is determined by the relative C configura- 
tion between coupled vortices, which energy is itself determined by the 
C-dependent direction of the stray fields induced by vortex-core shifts. 

Energy exchange rate and its dependence on interdistance. The 

energy exchange rate T ex between oscillators is technologically 
important from the information- signal transport point of view. It 
is determined by the coupling strength. Let T ex be the time period 
required for transferring the potential energy stored in disk 1 (here, 
due to the initial vortex- core shift) completely to disk 2. The value of 
T ex can be estimated from the frequency splitting determined from 
the modulation envelopes of |X X | and |X 2 |, (see Supplementary F) 
and consequently, is given as half of the modulation envelope period, 
^(27T/AcL>),sothatTex = 26 ± 3 nsand45 ± 6 ns for d int /(2R) = 1.05 
and 1.10, respectively. Since the strength of the interaction between 
two vortices determines the frequency splitting of dipolar-coupled 
vortices, a smaller value of T ex can be obtained by increasing the 
strength of the dipolar interaction. For any given material and 
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vortex state, it is known that the dipolar interaction between two 
vortices depends strongly on the interdistance between them 1415 . 
This was confirmed by the present experiments. Also, the 
micromagnetic simulation data (see Supplementary G) clearly 
show that the value of T ex varies with the interdistance as x ex 
~(d int /2R)" 3 91±007 . Accordingly, the energy transfer rate between 
vortex- state disks is tunable varying their interdistance. 

Energy exchange rate depending on relative vortex polarization 
configuration. It is worth noting that the relative motion of the 
vortex- core gyration of two vortices also affects the strength of 
their interactions: the dipolar interaction of vortex- core gyrations 
having opposite rotational senses due to opposite polarities is 
much stronger than that of gyrations having the same rotational 
sense due to the same polarity. From the simulation (see 
Supplementary H), T ex for the opposite polarity is ~ 2 times faster 
than that for the same polarity. Since the magnetostatic interaction 
strength depends directly on the stray field of each disk, thus, the 
larger the saturation magnetization is, the stronger will be the 
interaction strength, yielding faster T ex and hence faster signal 
transport. Because the polarity of a given vortex can be 
manipulated in a controllable manner by applying rotating or 
pulse fields locally to that vortex state using specially designed 
electrodes, as reported in Refs. [30-31], the polarities of vortices 
can be one of the crucial parameters for governing the interaction 
strength of energy transfer, and hence the exchange rate. 

Discussion 

To conclude, we experimentally observed, by time- resolved MTXM, 
energy and information-signal transfer via stimulated vortex gyra- 
tion through dipolar interaction between separated magnetic disks. 
This robust new mechanism for energy transfer provides the advan- 
tages of a fast and tunable energy transfer rate that is a function of 
disk interdistance and interaction strength. Control of energy loss 
during gyration-mediated signal transfer is possible with material 
engineering; in fact, almost lossless energy transfer can be achieved 
by employment of a material having negligible damping. Vortex 
gyration also can be achieved with low-power consumption through 
the resonant vortex excitation. This finding opens a new avenue to 
the development of energy efficient information processing devices 
such as logic gates based on vortex-state networks. 

Methods 

Sample preparation. The sample shown in Fig. 1 was prepared on a 100 nm-thick 
silicon nitride membrane by electron-beam lithography, thermal evaporation, and 
lift-off processing. The two disks in each pair have a different center-to-center 
interdistance, such that d int /2R = 1.05,1.1,1.15, and 1 .2, and each pair is separated at a 
sufficiently large distance, 5.6 /mi, from neighbouring pairs. A single-strip Cu 
electrode leads to strong local fields, as shown in Supplementary I, so that local 
excitations of the vortex, in each pair, are possible only in disk 1. Note that we 
confirmed that such local fields do not allow for vortex excitations in disk 2 positioned 
sufficiently far from the electrode stripline, in any of the d int /2R cases studied here. 
The eigenfrequency of the vortex gyrations in isolated Py disks, by ferromagnetic 
resonance measurements performed on an array of Py disk pairs of the same 
dimensions and d int I 2R = 1.56, was determined to be around 157 ± 3 MHz. 

Micromagnetic simulation. We performed micromagnetic simulations of coupled 
vortex gyrations under free relaxation for two Py disks of the same dimensions as 
those of the sample for d int /(2R) = 1.05 and 1.10. To mimic the experimental 
conditions, we used the same initial vortex-state configurations as those of the 
real sample. The material parameters for Py are as follows: the exchange stiffness 
A ex = 13 pj/m, the saturation magnetization M S =7.2X 10 5 A/m, and a zero magnetic 
anisotropy constant. The cell size was 4X4X50 nm 3 with the Gilbert damping 
constant a = 0.01. We used the OOMMF code that utilizes the Landau-Lifshitz- 
Gilbert equation of motion 32 . 
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